TUM-HEP-458/02 
hep-ph/0204225 
April 2002 



Enhancement of B (B d — ► ^i + ^-)/B (B s — ► in the MSSM 

with Modified Minimal Flavour Violation and Large tan/3 



C. BOBETH, T. EWERTH, F. KRUGER, AND J. URBAN* 
Physik Department, Technische Universitat Miinchen, D-85748 Garching, Germany 



Abstract 

We extend our previous analysis of the decays Bd,s i n the Minimal 

Supersymmetric Standard Model (MSSM) to include gluino and neutralino 
contributions. We provide analytic formulae, valid at large values of tan (3, 
for the scalar and pseudoscalar Wilson coefficients arising from neutral Higgs 
boson exchange diagrams with gluinos and neutralinos. Together with the re- 
maining contributions (VF ± , H^, x 1 * 1 ), and assuming the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix to be the only source of flavour violation, we assess 
their implications for the branching fractions B (Bd,s ^ + 1 1 )- Of particular 
interest is the quantity R = B (Bd - > fJ- )/S (B s — > since (i) the 

theoretical errors cancel to a large extent, and (ii) it offers a theoretically clean 
way of extracting the ratio |Vtd/Vt s | in the Standard Model, which predicts 
Rsm ~ \V t d/V ts \ 2 ~ O(10" 2 ). Exploring three different scenarios of modified 
minimal flavour violation (MFV), we find that part of the MSSM parameter 
space can accommodate large Bd,s branching fractions, while being 

consistent with various experimental constraints. More importantly, we show 
that the ratio R can be as large as O(l), while the individual branching frac- 
tions may be amenable to detection by ongoing experiments. We conclude 
that within the MSSM with large tan/3 the decay rates of Bd — > /i + ^ and 
B s — > can be of comparable size even in the case where flavour violation 

is due solely to the CKM matrix. 
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I. INTRODUCTION 



In the Standard Model (SM), the dominant contributions to the decays B q — ► / + /~, 
where q = d,s and I — e,fi, r, arise from Z°-penguin diagrams and box diagrams involving 
W ± bosons. The contributions due to neutral Higgs boson exchange, on the other hand, are 
suppressed by a factor of m^m^/M^ < 10~ 3 , and therefore are completely negligible. 

The SM predicts the B s — > ji + n~ branching ratio to be |],|| 

B{B S -> //+/0 = (3.2 ± 1.5) x KT 9 , (1.1) 
and the ratio of branching fractions 

^ r Bd M B J 2 Bd \V td \ 2 ^. 2 , 



B (B s -> 



SM *b. M Bs /i s \V U \ 2 



O(IO^), (1.2) 



where tr is the lifetime of the I? g meson, Mb ? and are the corresponding mass and 
decay constant. We note that the main uncertainty on the branching ratio in Eq. ( |1 . 1| ) arises 
from fs s , with a typical value of (230 ± 30) MeV obtained from lattice QCD calculations 
0. On the other hand, the relative rates of Bd and B s decays [Eq. ( |1.2|) 1 have a smaller 
theoretical uncertainty due to the appearance of the ratio f'B d /fB s , which can be more 
precisely determined than fg a alone. A determination of R can therefore potentially provide 
valuable information on |V^/V^| in the SM. However, given the SM prediction of BiB^ 
/i + /i _ ) ~ O(10 -10 ), the &d — > decay is experimentally remote unless it is significantly 

enhanced by new physics. Thus, the purely leptonic decays of neutral B mesons provide 
an ideal testing ground for physics outside the SM, with the current experimental upper 
bounds H 

B {B d -> (j.+fj,-) < 2.8 x 10~ 7 (90% C.L.), (1.3a) 



B (B s -> ti+fT) < 2.0 x 10~ b (90% C.L.). (1.3b) 

Our main interest in this paper is in a qualitative comparison of the Bd. s — * 
branching fractions in the presence of non-standard interactions, which can be made by 
using the ratio 

B {B s n+fM ) 

Referring to Eq. ( |1.2j ), it is important to note that the suppression of R in the SM is largely 
due to the ratio of the Cabibbo-Kobayashi-Maskawa (CKM) elements. This dependence on 
the CKM factors allegedly pertains to all models in which the quark mixing matrix is the 
only source of flavour violation. It is therefore interesting to ask if R could be of the order 
unity in some non-standard models where flavour violation is governed exclusively by the 
CKM matrix. Working in the framework of the Minimal Supersymmetric Standard Model 
(MSSM) with a large ratio of Higgs vacuum expectation values, tan f3 (ranging from 40 to 
60), we show that such a scenario does exist, and study its consequences for the Bd,s 
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branching ratios|] To this end, we extend our previous analysis to include the effects 
of gluino and neutralino contributions, in addition to those coming from diagrams with 
W ± , x ± , charged and neutral Higgs boson exchange. 

The outline of this paper is as follows. In Sec. [0], we define modified minimal flavour 
violation (MFV) and discuss briefly three distinct scenarios within the MSSM. The effective 
Hamiltonian describing the decays Bd )S i n the presence of non-SM interactions 

is given in Sec. |TJ, together with the branching ratios. In Sec. [TV], we present analytic 
formulae for the gluino and neutralino contributions to scalar and pseudoscalar Wilson 
coefficients governing the h — ► ql + l~ transition. Numerical results for the branching fractions 
£> (Bd )S an d the ratio R are presented in Sec. |V[ In Sec. |VT|, we summarize and 

conclude. 



II. THE MSSM WITH MODIFIED MINIMAL FLAVOUR VIOLATION 



There exists no unique definition of minimal flavour violation (MFV) in the literature 
(see, e.g., Refs. [P| P~2HT5|] ) . The common feature of these MFV definitions is that flavour 
violation and/or flavour-changing neutral current (FCNC) processes are entirely governed 
by the CKM matrix. On the other hand, they differ, for example, by the following additional 
assumptions: 



i) there are no new operators present, in addition to those of the SM JT|,||,|T2[|T3] . 

ii) FCNC processes are proportional to the same combination of CKM elements as in the 



SM [0 



iii) flavour transitions occur only in charged currents at tree level |L5 



While these ad hoc assumptions are useful for certain considerations, such as the construction 
of the universal unitarity triangle fl3f , they cannot be justified by symmetry arguments on 
the level of the Lagrangian. For example, the number of operators with a certain dimension 
is always fixed by the symmetry of the low-energy effective theory. Whether the Wilson 
coefficients are negligible or not, depends crucially on the model considered and on the part 
of the parameter space. Furthermore, the requirement that FCNC processes are proportional 
to the same combination of CKM elements as in the SM fails, for example, in the MSSM 
and can be retained only after further simplifying assumptions. The last statement in iii) is 
of pure phenomenological relevance in order to avoid huge contributions to FCNC processes. 

Using symmetry arguments, we propose an approach that relies only on the key ingre- 
dient of the MFV definitions in Refs. & 12-15 , without considering the above mentioned 



additional assumptions i) — iii) - We call an extension of the SM a modified minimal flavour- 
violating (MFV) model if and only if FCNC processes or flavour violation are entirely ruled 



1 A number of authors J|, ||-fj]| have noted that in the large tan/3 regime neutral Higgs-boson 
contributions can increase the SM decay rates of Bd,s —* by up to several orders of magnitude. 

Moreover, there is a correlation between B (B s — > /x + /i _ ) and Bg~B® mixing MM, and, in mSUGRA 
scenarios, between B (B s — > fi + fi~) and (g 



2) , |1(J (see also Ref. 
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by the CKM matrix; that is, we require that FCNC processes vanish to all orders in pertur- 
bation theory in the limit Vqkm - > 1- 

Let us briefly motivate our definition of MFV. The structure of the SM is such that the 
conservation of lepton numbers L e , and L T is exact. This is reflected by three global U(l) 
symmetries (one for each lepton flavour) of the SM Lagrangian. As a consequence, there are 
no FCNC processes in the lepton sector to all orders in perturbation theory. We regard the 
breaking of these lepton flavour U(l) symmetries as flavour violation.^ As far as the quark 
sector of the SM is concerned, the CKM matrix explicitly breaks the quark flavour U(l) 
symmetries. Thus, in the quark sector, FCNC processes can be generated in electroweak 
interactions. In the limit Vqkm ~^ 1 the three U(l) quark symmetries are restored and 
therefore FCNC transitions vanish to all orders in perturbation theory. This observation 
leads us to the definition of MFV given above. As will become clear, the advantage of MFV 
is that it is less restrictive than MFV, while the CKM matrix remains the only source of 
FCNC transitions. 

Our definition of MFV is manifest basis independent. However, in order to find a useful 
classification of different MFV scenarios within the MSSM, we will work in the super-CKM 
basis (see Ref. |17]] for details). In this basis the quark mass matrices are diagonal, and both 



quarks and squarks are rotated simultaneously. The scalar quark mass-squared matrices in 
this basis have the structure 

2 _ (Ml Ml\ 2 _ (Ml M 2 D \ 

u ~ I M 2] M 2 r D ~ I M 2] M 2 I ' ^ ' 

V JVI u LR JVI u RR J V JVI d lr JVI d rr J 

where the 3x3 submatrices are given by 

M 2 Ull = Ml L + My + -M 2 Z cos 2/3(3 - 4 sin 2 9 w )t, (2.2a) 

M 2 Ulr = MuiAh - //cot/311), (2.2b) 
M 2 Urr = Ml R + M 2 + 2 -M 2 z cos 2/3 sin 2 6 W H, (2.2c) 



M 2 D =M\ +M 2 D - -Ml cos 2/3(3 -2 sin 2 6 W ) 11, (2.3a) 



M 2 Dlr = M D (A* D - /t tan/31), (2.3b) 



M d rr = M 2 6r + M 2 D - \m 2 z cos 2/3 sin 2 9 W TL. (2.3c) 



2 This is different from Ref. [16], where flavour violation is regarded as the breaking of the flavour 
symmetry Gf- 
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Here, M~ and M? are the soft SUSY breaking squark mass-squared matrices, \i is 
the Higgsino mixing parameter, Ay and A u are soft SUSY breaking trilinear couplings, 1 
denotes the 3x3 unit matrix, and 

M v = diag(m u , m c , m t ), M D = diag(m d , m s , m h ). (2.4) 

Because of SU(2) gauge invariance, the mass matrix is intimately connected to M?. 
via 

Ml =Vi KM MlV CKM , (2.5) 



which is important for our subsequent discussion. In the case of MFV, the matrices Ay, A u , 
M% , M~ must be diagonal: 

Ur Ur ^ 

A v = diag(A„, A c , A t ), A D = diag(A d , A„ A b ), (2.6) 



M l R = diag(ml fl , mj R , m|J, M\ r ee diag^, m\ R , m?J. (2.7) 

Furthermore, assuming that there are no new CP-violating phases, in addition to the single 
CKM phase, the matrices Ay and A u , as well as /i, are real. 

Taking into account the relation in Eq. (|2.5|), one encounters three cases of MFV. 



Scenario (A): 

M v is proportional to the unit matrix, and so M?. = . As a result, there are no 
gluino and neutralino contributions to flavour-changing transitions at one-loop level. 
The implications of this scenario for rare B decays have been discussed, e.g., in Refs. || 
^-0. This scenario of MFV coincides with the MFV scenario at low tan/?, as defined 
in Refs. OII1. 



Scenario (B): 

is diagonal but not proportional to the unit matrix and, in consequence, M?. 
has non-diagonal entries. In such a case, there are again no gluino and neutralino 
contributions to flavour-changing one-loop transitions involving only external down- 
type quarks and leptons. However, additional chargino contributions show up, due to 
non-diagonal entries of M?. . This scenario has recently been investigated in Ref. 0. 



Scenario (C): 

M? is diagona 

diagonal entries in . Accordingly, gluino and neutralino exchange diagrams (in 

lvl 

flavour- changing transitions at one-loop level that involve external down-type quarks. 



M?. is diagonal but not proportional to the unit matrix, which gives rise to off- 
anal entries in M\ 

addition to those involving W ,x , charged and neutral Higgs bosons) contribute to 



The common feature of all these scenarios is that the CKM matrix is the only source of 
flavour violation. For the remainder of this paper we will concentrate mainly on scenarios 
(B) and (C). To our knowledge, the consequences of scenario (C) for the decays B d s — > ji + /i~ 
have not yet been discussed in the literature. 
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III. EFFECTIVE HAMILTONIAN AND BRANCHING RATIO 



The effective Hamiltonian responsible for the processes B q — > / + / , with q = d,s and 
I = e, fi, t, in the presence of non-standard interactions is given byf\ 

Hex = -^V tb V t * q £ [ci^Oi^ + cJ^QKA*)], (3-1) 

with the short-distance coefficients c['\fi) and the local operators 

do = (q^PMh^l), O' 10 = (qrP R b)(h^ 5 l), (3.2a) 

O s = m b (qP R b)(ll), 0' s = m q (qP L b)(ll), (3.2b) 

Op = m b (qP R b)(h 5 l), 0' P = m q (qP L b)(h 5 l), (3.2c) 

where Pl,r — (1 T 7s)/2- In addition to the operators in Eqs. (|3T2"|), there are tensor 
and vector operators, (qa^ Pl^X/ct^P^^) and {q r ) fl Pi J fib){l r )^l). However, they do not 
contribute to the B q — > decays and do not mix with those appearing in Eq. Q3.ip . 
As a matter of fact, the matrix element that involves the antisymmetric tensor o^ u must 
vanish since p^ = p B is the only four-momentum vector available. Similarly, the matrix 
element (/ + /~ 1/7^/10) does not contribute when contracted with {Q\q^P L ^ R b\B q {p)) cx p M by 
the equation of motion. 

Evolution of the short- distance coefficients c['\fi) from the matching scale /i 4 = m\ ole 
down to the low-energy scale fib = m b° le can t> e performed by means of the renormalization 
group equation. We note that the anomalous dimensions of the operators in Eqs. ( p.2| ) 
vanish, and thus the renormalization group evolution is trivial. 

Because of the pseudoscalar nature of the decaying B q meson, the hadronic matrix ele- 
ments are non-zero only for axial- vector and pseudoscalar operators; namely, 

(0|97/i75&|-B 9 (p)) = ip»fB q , (3.3) 
and, employing the equation of motion, 

(0\q l5 b\B q (p)) = -if Bq ^ ■ ( 3 - 4 ) 

q m b + m q 

Here f Bq is the B q meson decay constant, which can be obtained from lattice QCD compu- 
tations ||: 

f Bd = (200 ± 30) MeV, f Bs = (230 ± 30) MeV, J -p- = 1.16 ± 0.04. (3.5) 

JB d 



3 We omit the index q on the operators and the corresponding short-distance coefficients. 
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(Similar results have recently been obtained from a QCD sum rule analysis fl8|| .) 

If the lepton spins are not measured, the branching ratio for the case I = fi takes the 
general form 

GWMl f 2 t b 



64tt 3 



v th v, 



* 1 2 



tb v tq\ 



1 - Ami 

(J, 



l-Afn z u )\F s \ z + \F P + 2m^F M 



with the notation = m^/MB q and the dimensionless form factors 



S,P 



M, 



B„ 



cs,pm b 



z s,P m q 



ClO 



'10- 



(3.6) 



(3.7) 



m b + m q 

In the SM, the contributions involving the neutral Higgs boson are completely negligible, 
and so B (B q — > /U + /i~) oc m 2 , which is a consequence of helicity suppression .[J] Finally, 
allowing the input parameters in Eq. (|3.6|) to vary over the interval in Eq. ( |3.5|) , and using 
the ranges for the CKM factors given in Ref. jl], the ratio of decay rates of Bd,s ^ + fi~ 
within the SM is expected to be in the range 



0.02 < R SM < 0.05, 
which is largely due to the imprecisely known CKM elements. 



(3.8) 



IV. HIGGS-BOSON CONTRIBUTIONS TO THE DECAYS B d , s -> l+l~ 

We now turn to the computation of the scalar and pseudoscalar Wilson coefficients in 
the b — > ql + l~ transition arising from gluino and neutralino exchange diagrams within the 
general MSSM. As mentioned earlier, we perform our calculation in the large tan/? regime 
(i.e. 40 ^ tan/3 ^ 60). For the remaining contributions (W ± ,H ± ,x ± ) to these short- 
distance coefficients, we refer to Refs. @, • 

The relevant box and penguin diagrams are depicted in Fig. [I], where H°, h°, A and 
G° are the neutral Higgs and would-be-Goldstone bosons respectively, l a are the charged 
sleptons, d a denote the down-type squarks, x°k are the neutralinos, and g represents the 
gluino. We perform the calculation in the 't Hooft-Feynman gauge, using the Feynman rules 
of Ref. |19[ , and adopting the on-shell renormalization prescription described in Ref. . 



In our subsequent calculation, we will exploit the tree-level relations 

M% =M 2 H - M 2 H0 + M 2 h0 = M% + Mf, (4.1) 

sin2a _ ( M 2 H , +M fe 2 \ cos 2a _ / M\ - M 2 \ 

sin2/? = \M 2 H0 - M 2 h J ' cos 2(3 = \M 2 H0 -M 2 h J' { ' ' 

where M40 and Mh are the masses of the CP-odd and charged Higgs boson respectively. 
M h o H o and a are the masses and mixing angle in the CP-even Higgs sector. This leaves two 
free parameters in the Higgs sector that we choose to be Mh and tan j3. 



Because the B meson has spin zero, both fj, + and fi must have the same helicity. 
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d, s 



d, s 



d, s 



xV 



Ixl 



xl 



xl 



xl 



to uo m no 



d, s 



d a \ / d b 

V 



I ,h ,A°,G 



d, s 



V 



H°,h°,A ,G 



d, s 




FIG. 1. Gluino and neutralino diagrams contributing to the decay B^^ — > l + l~ , with I = e,fx, t. 
The indices k, k' and a, b run from 1 to 4 and from 1 to 6 respectively. The Feynman rules for the 
counterterm diagrams (last diagram) can be found in Eqs. (5.1) and (5.2) of Ref. M. 



Introducing, for convenience, the dimensionless variables 

,2 



2 

m- 

~ M 2 "' Xak 



m 



d a 



2 ' 



(4.3) 



our results for the gluino and neutralino contributions at large tan f3 can be summarized as 
follows. 

(i) Gluino: 



1 4a. 



mi tan 2 f3 



^S,P 



V tb V t * 3a m b (M 2 H - M^)M~ 9 



£(r^) ?6 (r^) 



«3 



a,b=l 



X 



D 2 (x a , x b ) [T D ^M D A* D T D ^ ± H.c.] 6a - M? g D 3 (x a )5 ab ] 



(4.4) 



'Our results for the gluino and neutralino contributions differ somewhat from those given in 



Ref. [20 1, which contains typographical errors in the formulae for the various Wilson coefficients [21]. 
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>~ 9 _ 1 4a s m,tan 2 /3 ^ Dr] d 
S ' P ~ V tb V* 3a m q {Ml - M^)M- g 2-^ W )a3 

x {£> 2 (x ,x b ) [T Dl M d A* d T d ^ ± H.c.] 6a =F Mf D 3 (za)<u}, (4.5) 

where the subscript q on the T matrices (see Appendix [FJ for details) is equal to 1 (2) 
for B d (B s ) decays. Furthermore, 9 W is the Weinberg angle and the functions -D 2 ,3 are 
given in Appendix Observe that in Eqs. ( |4.4| ) and ( |4.5| ) we have retained the leading and 
subleading terms in tan (3, which might become comparable in size in some part of the MSSM 
parameter space. A term is called leading or subleading if it has the following structure: 

^| leading ~tan^/3(-) , ^ suWeading ~tan^ (-) , (4.6) 

n — 0, 1, 2,. . . . Further, m denotes lepton and light quark masses while M stands for masses 
of particles that have been integrated out. 
fii) Neutralino: 



° S ' P = V tb V t ll2M w {Ml-Ml)^9 w h h %l 



V tb V t *12M w (M 2 H - M^)sin 2 9 
x { tan 4 ^NU^UT^U 



3(Xak 



k=l a=l 



+ tan 3 /3N k3 



2?n 

(tan9 w N kl - 3N k2 )(T D ^) qa {T D ^) a3 + tan 9 w N kl (T D ^) qa (T D R ) a3 



+ tan 2 /5^tan^iV fcl {tnn9 w N kl - W k2 ) (T D ^) qa (T DR ) a3 ) , 
3m b J 



(4.7) 



mi 



x 



V tb V tq Y1M W {MI - M 2 ,) sin 2 9 W 



tan 4 P^NS^%^U 



4 6, 
fc=l a=l 



+ tan 3 [3N k * 3 



2mi 



(tan^iV^ - 3Nl 2 )(T D ^) qa (T D ^) a3 + ^tan^A^r^)^* 



'ci3 



+ tan 2 /5^tan^iV* 1 (tan^A^ - 3iV fc * 2 ) (r D ^ )(?a (r^; 



a-3 



(4.8) 



where AT is the neutralino mixing matrix, defined in Appendix [B]. Unlike the gluino con- 
tributions, we have kept only the leading term in tan/? [cf. Eq. ( |4.6| )1, which stems from 
the counterterm of the electroweak wave function renormalization (see Fig. [l]). However, we 
have checked numerically whether the subleading term is important. 
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V. NUMERICAL ANALYSIS 



A. Experimental constraints 



Currently available data on rare B decays already constrain new-physics contributions 
to the various Wilson coefficients governing the h — > s transitions.^ We take into account 



the following experimental upper limits [23,24 



B (B -> X s e + e~) < 10.1 x 10" 6 (90% C.L.), 



(5.1) 



B (B -> X s /i + ii-) < 19.1 x 10~ 6 (90% C.L.), 



(5.2) 



B {B -> K*n + fi-) < 3.1 x 10~ 6 (90% C.L.), (5.3) 
besides those on B d s — > /i + /i~ decays given in Eqs. (|1.3|). As far as the measured branching 



fractions of B — > Kfi + /j, [24] and B — > X s 7 [25] are concerned, we will allow the following 
ranges: 



0.5 x 10~ 6 ^B(B K/j+iT) sC 1.5 x 10" 



(5.4) 



2.0 x 10" 4 ^ B (5 -> X s7 ) ^ 4.5 x 10~ 4 . 



(5.5) 



The constraints on the flavour-changing entries in the matrices and will be dis- 
cussed below. 



B. Implications of MFV for the decays Bd, s — ► A 4 A 4 
In our analysis, we will perform a scan over the following ranges of MSSM parameters: 

125 GeV M H ^ 500 GeV, (5.6a) 



100 GeV ^ Mi 2 ^ 1 TeV, 250 GeV < M~ g < 1 TeV, 



(5.6b) 



1 TeV ^ {n, (Au)u, (A D )ii} ^ 1 TeV, 



(5.6c) 



300 GeV ^ {(Mjy L )ii, (M 3 J ih (M &R ) ii: (M bR ) tt ) <: 1 TeV. 



(5.6d) 



6 A recent analysis of exclusive and inclusive rare B decays within SUSY has been performed in 
Ref. [22] by using the SM operator basis. Recall that in the present work we utilize an extended 
operator basis. 
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Further, we fix tan/3 = 50 and m ; ~~ = 100 GeV, and take the lower bounds on the sparticle 
masses from As far as the decay constants fg are concerned, we will take the central 
values as given in Eq. ( |3.5| ). The remaining input parameters are fixed as given in Tables 1 
and 2 of Ref. 0. 

We have analyzed more than 100,000 parameter sets, which have been produced ran- 
domly. The so-called unconstrained plots contain those parameter sets that are consistent 
with the lower bounds on the sparticle masses. On the other hand, in the constrained plots 
we have also taken into account the bounds from Eqs. (|5.1| )- (|5.5| ) as well as the constraints 



on the flavour-changing entries in the matrices and (see discussion below). 

As for the various branching ratios, it turns out that the primed Wilson coefficients are 
negligibly small, and hence can be safely neglected. [Comments on this issue can be found 
at the end of the discussion of scenario (C).] 



Scenario (A) 



Recall that in scenario (A) the matrices M?. and are equal and proportional to 
the unit matrix. Therefore, the gluinos and neutralinos do not contribute at one-loop level. 
In Fig. |2], we have plotted B (Bd — > versus B (B s — > where the dashed line 

represents our reference curve with the slope Rsm while the solid line is our result obtained 
within scenario (A). The scan over the parameter region given in Eqs. Q5.6| ) shows that the 
ratio R is approximately constant and close to Rsm ~ 0.03. 




0.5 1 1.5 2 

B(B S -> n+fT) x 10 6 



FIG. 2. Predictions for the branching ratios B (Bd — * fj, + fj, ) vs. B (B s — > ) in scenario 

(A). The dashed line corresponds to Rsm ~ 0.03 while the solid line is the result for this scenario. 
The ratio R = B (B d -> fi + fi-)/B(B s -» varies between 0.028 < R < 0.029. Since the 

results for the constrained and unconstrained plots are indistinguishable, we only show the former. 
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Scenario (B) 



In this scenario, the matrix is diagonal, 



diag(m| m^mf 



(5.7) 



with at least two different entries; hence there are no gluino and neutralino contributions. 
Employing the relation in Eq. (|2.5| ), the matrix M?. becomes 



M 2 0l = VckmMIV^ 



/ m 2 - A 1 - 2 A 1 - 3 \ 

12* 2 A23 
Ul '"- ■ Al 



A 13* A 23* 2 



(5.8) 



with the off-diagonal elements 



A£ = (m 2 - m?)V^£ + (ml - m 2 s )V ub V; b , 



u L 



(5.9a) 



A 1 - 3 

u L 



(ml -mjWjVZ + im 



j 

6 £ 



(5.9b) 



)v cd v; d + 



m 



SL 



)v cb v*. 



(5.9c) 



In writing these equations, we have used the unitarity of the CKM matrix. The flavour- 
changing entries given above, together with the corresponding elements in the down squark 
sector [see Eqs. ( |5.12[ ) below], are constrained by experimental data on K°-K°, B°-B°, 
D°-D° oscillations, and the b -> 57 decay 0,|2|,|9|.[] 

In Fig. [| we have plotted B (B d — ► /i + /i _ ) versus B (B s — > /i + /U~) for unconstrained and 
constrained parameter sets. In case of the former, the deviation of R from R^m can be one 
order of magnitude, with the range 0.002 < R < 0.115. Applying the constraints on the 
flavour- changing entries in the matrix , as well as the bounds given in Eqs. fl5.1| )-( |5"3D , 
our numerical predictions for R in scenario (B) are similar to those in scenario (A). In fact, R 
varies between 0.026 and 0.030. It is important to note that the bounds on A*-? |[7J severely 
constrain the additional chargino contributions in scenario (B) in contrast to scenario (A). 



C. Scenario (C) 

We now repeat the analysis of the Bd yS — > [i + fi~ decays within the framework of scenario 
(C). In this case, the matrix is diagonal, 



M? l = diag(m| r , m\ , m\ ) 



t L ' 



(5.10) 



7 Strictly speaking, these constraints are only valid if all squark masses are close in size, but the 



order of magnitude should also be valid for non-degenerate masses [ 30 1 
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FIG. 3. Predictions for the branching ratios B (Bd — > /x + ^~) vs. B (B s — > /U + /i~) in 
scenario (B). The left plot shows the unconstrained parameter sets, with varying between 
0.002 < i? < 0.115. The right plot exhibits the results for the branching fractions, using the 
constrained parameter sets. In this case, R varies between 0.026 and 0.030. 



with at least two different entries. According to the relation in Eq. 



, this implies that 



M? has non-diagonal entries, so that gluinos and neutralinos contribute to the b — > ql + l~ 
transition already at one-loop level. In this case, can be written as 



A 12* m 2 A 23 



13* 



L 

23* 
Dl 



m 



2 



(5.11) 



where the flavour- changing off-diagonal entries are given by 



A 



(ml - ml T )V us V* d + {ml - m\ T )V ts V t 



tdi 



(5.12a) 



A 1 , 3 = (ml L - ml L )V ub V: d + (m? - m\ h )V th V t 



Dr 



[171; 



tdi 



(5.12b) 



(mj t - m c 2 r )K fc C + (mf - m?J^. 



(5.12c) 



As before, we take the constraints of Refs. |T7|,|28|,f29| on these off-diagonal elements 
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0.5 1 10 100 0.5 1 10 100 

B{B S -+ n+fi-) x 10 s B{B t -+ /i+ju-) x 10 s 



FIG. 4. Predictions for the branching ratios B (Bd — ► /U + /i~) vs. S (S s — > [i + (upper plots) 
and i? vs. (B s — > [i + (lower plots) in scenario (C). The left plots show the unconstrained 
parameter sets, with R varying between 0.026 < R < 20.636. The right plots exhibit the results 
for the branching fractions and their ratio R, using the constrained parameter sets. In this case, 
R varies between 0.026 and 2.863. 
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The scatter plots in Fig. f| exhibit an order-of-magnitude deviation from Rsm ~ 0.03. In 
the unconstrained case (left plots) R is predicted to be in the range 0.026 < R < 20.636, 
while for the constrained parameter sets (right plots) we find 0.026 < R < 2.863. A 
noticeable feature of scenario (C) is that there exists a lower bound on R, i.e. R > 0.95 Rsm 
(see upper plots of Fig. [|), which is due to the structure of the CKM matrix. We stress that 
this bound is valid only within scenario (C) and does not apply to scenario (B) or scenarios 
with new sources of flavour violation (see Sec. |TJ). 
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FIG. 5. B{Bd,s / U+ A t_ ) as function of the charged Higgs boson mass, Mh, for tan/3 = 50 
(left plot) and tan/3 = 60 (right plot), taking into account the experimental constraints on the rare 
B decays and on the flavour-changing entries. The remaining parameters have been kept fixed, as 
described in the text. 



In Fig. |5], the branching ratios of — » /i + yT and B s — > fi + n~ are displayed as functions 
of the charged Higgs boson mass, Mh, for tan/3 equal to 50 (left plot) and 60 (right plot). 
The remaining parameters are given by 



/x = -800 GeV, Mt = 500 GeV, M 2 = 200 GeV, M~ g = 250 GeV, 



(5.13a) 



A U = A D = diag(100 GeV, 100 GeV, 75 GeV), 



(5.13b) 



M u R 



M D R 



diag(500 GeV, 500 GeV, 500 GeV), 



(5.13c) 



M u L 



diag(460 GeV, 500 GeV, 517 GeV). 



(5.13d) 



Both plots are consistent with the constraints on rare B decays in Eqs. (|5.1| )- ([5T5| ), and 
with those on the flavour- changing entries []17| . Interestingly, in the left (right) plot, 
the ratio R ranges between 0.15 < R < 0.81 (0.44 < R < 1.78), while the magnitude of 
the individual branching fractions decreases drastically with increasing charged Higgs boson 
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mass, Mh- Note that for small Mh and tan/? close to 60 both branching ratios are in a 
region that can be probed experimentally, in Run II of the Fermilab Tevatron, BABAR and 
Belle. The large enhancement of R is due to Im [c? sp (Bd — > being by an order-of- 

magnitude larger than Im [c g SP (B s — > /i + |U~)] and a cancellation between the real parts of 

c| P and c§ P . As a consequence, the neutralino contributions to the scalar and pseudoscalar 
coefficients, although negligible compared to the corresponding contributions in the chargino 
and gluino sector, cannot be neglected in the B s — > u, + ur decay. The neutralinos contribute 
significantly to the ratio R; for illustration, the ratio R ~ 0.94 in the right plot (tan/3 = 60) 
at Mh = 130 GeV increases to R ~ 1.55 after neglecting the contributions from neutralinos. 

As already mentioned, we have found that the primed Wilson coefficients, 
numerically suppressed. In fact, applying the mass insertion approximation as outlined in 
Appendix [C], one can convince oneself that the enhancement factor l/m q (recall q = s,d) 
appearing in the scalar and pseudoscalar coefficients cancels out. Therefore, the primed 
Wilson coefficients are suppressed by a factor m g /mb <C 1, compared to the unprimed 
coefficients. 

VI. SUMMARY AND CONCLUSIONS 

In the framework of the minimal supersymmetric extension of the SM with large tan /3, 
we have computed the gluino and neutralino exchange diagrams contributing to the purely 
leptonic decays Bd, s — ► At + /i~. Together with our previous analysis [Q], the present paper 
provides a complete study of the decays Bd,s — > fi + u,~ at one-loop level in the MSSM with 
modified minimal flavour violation and large tan (5. 

We have defined MFV using symmetry arguments and have shown that MFV is less 
restrictive than MFV, while the CKM matrix remains the only source of flavour violation. 
We have given a criterion for testing whether a theory belongs to the class of MFV. Within 
the MSSM we have investigated three scenarios that are possible within the context of MFV. 
In particular, we have studied the case where the gluino and neutralino exchange diagrams 
contribute besides W ± , H ± , X s1 [scenario (C)]. The neutralino Wilson coefficients are numer- 
ically smaller than those coming from the chargino and gluino contributions. However, we 
have found that in certain regions of the MSSM parameter space cancellations between the 
chargino and gluino coefficients occur, in which case the neutralino contributions become 
important. As a matter of fact, for the SUSY parameter sets examined, we found that a 
large value of R = B (B d — > n~)/B (B s — > always involves such a cancellation. 

Including current experimental data on rare B decays, as well as on K, B, D meson 
mixing, we found that in certain regions of the SUSY parameter space the branching ratios 
B (Bd — > and B (B s — ► can be up to the order of 10~ 7 and 10~ 6 respectively. 

Specifically, we showed that there exist regions in which the branching fractions of both 
decay modes are comparable in size, and may well be accessible to Run II of the Fermilab 
Tevatron. 

We wish to stress that a measurement of the branching ratios B(Bd, s / i+ A t ~) 5 or 
equivalently, a ratio R of 0(1), does not necessarily imply the existence of new flavour 
violation outside the CKM matrix. Nevertheless, any observation of these decay modes in 
ongoing and forthcoming experiments would be an unambiguous signal of new physics. 
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APPENDIX A: AUXILIARY FUNCTIONS 

The loop functions appearing in the formulae of Sec. [TV] are given by 

D 2 {x,y) = Xl * 1X + (x <-> y), (Al) 

(1 - x){x - y) 

D 3 (x) = - . (A2) 

I — x 



APPENDIX B: MASS MATRICES 
1. Squark mass matrices 

The squark mass-squared matrices given in Eq. ( |2.1| ) are diagonalized according to 



Y u MlY u ^ = diag^ 



T D M 2 D T D ^ = diag( 



di' 



(Bl) 



It is convenient to split the 6x6 squark mixing matrices, T U ' D , into two 6x3 submatrices: 

(r«) ai = (r^) ai , (r«) a , l+3 = (r^) oi , q = u,d, (B2) 



where a = 1, . . . , 6 and i = 1,2, 3. 



2. Neutralino mass matrix 



The neutralino mass matrix has the structure [[31 

/Mi -M z sin 6 W cos (3 M z sm9 w sm(3 \ 







Mo 







—M z sin 6 W cos (3 M z cos 9w cos (3 
\ M z sin 9w sin (3 —Mz cos 0w sin (3 —\i 

which is diagonalized by a unitary matrix iV such that 

N*M^N^ = diag(M^o,M^o,M^o,M^o 



M z cos 6w cos (3 —M z cos9w sin P 







(B3) 



(B4) 
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APPENDIX C: PERTURB ATIVE DIAGONALIZATION OF M 2 D 

Within the framework of scenario (C), all flavour-changing interactions and CP viola- 
tion are due exclusively to the CKM matrix. In this case, the matrices -M 2 Dlr and M-\> RR 
in Eq. (|2.1|) are diagonal and real whereas M- 2 d ll contains complex off-diagonal entries 
[cf. Eq. (|5.11| )1. In order to diagonalize the down squark mass-squared matrix M. 2 D pertur- 
batively, we rewrite it as 

Ml = Ml + A D= ( "h') + f V n\ (CD 



M'' Dlr MI r J V °. 

where M. 2 D (Ao LL ) contains the diagonal (off-diagonal) elements of Ai 2 D . We remark 
parenthetically that for large values of tan j3 the LR elements in Ai 2 D cannot in general be 
treated as small perturbations. [This is the reason for choosing the decomposition of A4 2 D 
given in Eq. (|CT|).] 

Writing T D as a product of two unitary matrices, T D = XY, we have 

(A^)diag = T D M 2 D Y D ^ = X {D + A D } X\ (C2) 

where D = YM 2 D Y^ = diag(£>i, . . . , As) and A D = YA D Yl Then, if |(A D )^| < D a , we 
can make the ansatz 

X = 11 + 5X {1) + 5X {2) + ■ • • , 
(^D)diag = D + SD® + SD^ + ■■■ (C3) 

to solve Eq. flC2] ) perturbatively in terms of Ad. As a result, we obtain 

T^/K-jrf, = Yif(D c )Y cb + YlA cd f x {D c1 D d )Y db 

+ YlA ce A ed f 2 (D c , D d , D e )Y db + • • • , (C4) 

where / denotes an arbitrary loop function, and 

tt \ fi x )-f(y) f t \ h{x,z) - fi(y,z) 

fi(x, y) = > h{x, V, z) = , (C5) 

x — y x — y 

for x ^ y A x ^ z A y ^ z. In all other cases the corresponding limit of the functions /i 2 has 
to be taken. Setting Y = 1L in Eq. (|U4|) reproduces the result given in Refs. 
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